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Abstract 

m 
> 

^^ • Considering the low-temperature (T) and high-baryon-number-density {ub) 

Qs ' region of the QCD phase diagram, we present a model for the first-order 

o' 

^^ ' phase transition between the quark-gluon plasma (QGP) and the recently 

O 

c~| ' proposed colour superconducting phase. We study nucleation of a droplet 

Q^' of the superconducting phase within the metastable QGP gas. Numerical 

1-^ ' results for the activation energy, radius and other physical parameters of the 

[^ • droplets, at various temperatures, densities and gap parameters, are given. 

X 

j^ ■ We have estimated the latent heat of the phase transition. In the T — ub 
plane, we are able to demarcate the region of the superconducting phase. 
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I. INTRODUCTION 

Our understanding of the QCD matter at low temperatures (T) and high baryon number 
densities (n^) has developed rapidly over the past few years; for recent reviews, see Jlj-Q. 
A likely scenario is that in this region, quarks form Cooper pairs and new condensates 
develop. One of the most interesting questions today concerns the detailed mapping of 
the QCD phase diagram, especially in the above region where QCD may exhibit a colour 
superconducting phase. In this context, the light quark masses {'mu,d,s) are expected to play 
an important role. Figure 1 shows two likely scenarios corresponding to rUs ^ Triu^d 7^ 
and rus = mu^d — 0. The transition between the quark-gluon plasma (QGP) phase and the 
colour superconducting phase is expected to be first order, in QCD |l|,@]. 

In this paper, we study the first-order phase transition between the QGP phase and the 
colour superconducting phase (Fig. 1). Our approach is phenomenological and is similar in 
spirit to that adopted by Cote and Kharchenko 0. They, however, studied the Bose-Einstein 
condensation in a gas of trapped atomic hydrogen at temperatures ~ 50 fiK and densities 
~ 10^^ cm~^. Temperatures and densities that we shall consider here are, of course, vastly 
different, and QGP has very little in common with the atomic hydrogen gas. Interestingly, 
however, as we shall see, the same physics of activation energy barrier governs the formation 
of droplets in both cases. 

It is important to get a handle on the cold and dense quark matter because of its possible 
relevance to the neutron-star core. While the region of asymptotic densities is amenable to 
rigorous calculations, the same cannot be said about the region of intermediate densities 
studied here. Both numerical (lattice) and direct experimental studies of this region of the 
QCD phase diagram are out of reach at the moment. Hence a phenomenological study such 
as the present one is appropriate and worthwhile. It has allowed us to make quantitative 
statements about the schematic phase diagram in Fig. 1. 

Alford et al. [0] have discussed droplet formation in the transition between the hadronic 
matter phase and the 2SC phase (Fig. 1), at T = 0, as a function of density. The chiral 



symmetry is restored inside their droplets (which they call nucleons) but is broken outside. 
Neergaard and Madsen ^ also have studied the droplet formation but in the QGP to 
hadronic matter transition occurring at low baryon number densities. In this paper, we 
have considered the QGP to colour-superconducting phase transition occurring at high ub- 
To our knowledge, droplet formation in this transition has not been discussed previously in 
the literature. Our droplets are the opposite of those in in the sense that the chiral 
symmetry is broken inside our droplet, but is restored outside. 

II. MODEL 

Consider a gas of weakly interacting quarks, antiquarks and gluons in a volume V. 
We take q and q to be the massless current quarks (Fig. 1). We are interested in low 
temperatures (T) and high number densities [n) such that nX^ ^ 1 where A is the thermal 
wave length (see Appendix). We also need the baryon number density ub = {uq — ng)/3 to 
be large compared to that for the normal nuclear matter, namely 0.17 fm~^. 

If the interactions among the particles are treated to the lowest order in a^, the strong 
coupling constant (Fig. 2), the baryon number density is given by (c = /i = 1) 
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where dq = 2 (spin) x 3 (colour) x 3 (flavour) = 18 is the quark degeneracy factor, k is the 
Boltzmann constant and fiq is the quark chemical potential 0. For a given baryon number 
density ub and temperature T, Eq. (1) uniquely fixes fiq. This in turn determines the total 
energy and pressure of the gas as 
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Here dg = 2 (spin) x 8 (colour) = 16 is the gluon degeneracy factor. It was possible to 
obtain the above analytic expressions because we have considered the gas of quarks as well 



as antiquarks, rather than a gas of quarks alone 0. At low T and high ns, the /ig terms 
in the pressure dominate and the gluonic pressure arising from the first term in Eq. (2), is 
quite negligible. 

We now consider the first-order phase transition from QGP to the superconducting phase 
as the temperature is lowered (Fig. 1). In the latter phase, as a result of the attractive one- 
gluon-exchange interaction between quarks in the colour 3 channel, quarks near the Fermi 
surface tend to pair up as bosons. These are the Cooper pairs which we assume to be 
massless. Density fluctuations in the metastable QGP gas may provide suitable conditions 
for nucleation and growth of a droplet of the superconducting phase. The high-temperature 
phase (QGP) is present around the droplet and the low-temperature superconducting phase 
exists inside the droplet. We will presently explore the energetics for the formation of such 
a droplet by calculating the energy barrier between the two phases. 

We have given above the expressions for number density, energy and pressure of the 
QGP; see Eqs. (l)-(3). We now derive corresponding expressions for the superconducting 
phase inside the droplet. Before we begin, it is necessary to recall important differences 
between the BCS superconductor and the colour superconductor, (a) In the BCS theory of 
superconductivity the ratio of the energy gap (A) to the Fermi kinetic energy for a typical 
metal is of the order of 10~"^. In contrast, in the problem under consideration, the gap A 
is estimated to vary between several tens of MeV to about 100 MeV |]l|-^, and the quark 
chemical potential fig appearing in Eq. (1) takes values between ~ 250 and ~ 600 MeV for 
the hb and T considered here, (b) In the BCS superconductor, the Debye energy uj£, is much 
smaller than the chemical potential /z, while in the colour superconductor, the two may be 
of the same order of magnitude, (c) Finally, if n is the density of electrons and ^ the Pippard 
coherence length which measures the spatial extension of the Cooper pair wave function, 
then in the BCS theory, n^'^ ^ 1. In the present case, on the other hand, ^ ~ l/(7rA) ~ 0.6 
fm for A = 100 MeV, and n^'^ < 1 for the densities considered in this paper. This suggests 
that the qq pairs are rather compact objects unlike Cooper pairs of electrons in a metal. For 
a good discussion of these and related issues, see a review article by Kerbikov [|l^ . We shall 



make use of the above points while deriving expressions for thermodynamics quantities of 
the droplet. For the sake of simplicity, we shall assume that all quarks in the energy interval 
(/ig — A, fiq) form Cooper pairs and all other quarks remain unpaired, and shall ignore other 
interactions. 

Consider a spherical region of radius Ra, volume V^, containing A^i quarks out of which N2 
remain unpaired and N3 undergo Cooper pairing, when the pairing interaction is "switched 
on" . We have 



N, = N,- -^ / 4V/^P, (5) 

iVs = iVi - N2, (6) 

where / = [exp /3(p — fig) + 1]~ with j3 = l/{kT). The corresponding energies -£"1,2,3 are 
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E2 = E,- -f^ / 4V/dp, (8) 

--(f)(i-)- (^) 

In the last equation {N^/2) is the number of Cooper pairs (bosons) in the droplet and {n/Rd) 
is the energy of the ground state of the single-particle spectrum of a spherical cavity of radius 
Rd- Due to the small size of the droplet, even the lowest excited state has energy which 
is very large compared to the temperature. Therefore all the qq pairs tend to occupy the 
ground state [Q. Next, entropies 5'i,2,3 are given by 
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^i = -^y^ Wdp[f\nf+{l-f)\n{l~f)], (10) 

52 = S, + ^r 4vrp>[/ln/ + (l-/)ln(l-/)], (11) 

53 = 0. (12) 

5*3 can be taken to be zero for the reasons given above. Finally, the pressures -Pi,2,3 are 
calculated as follows: 



Pi = Ej{Wd). (13) 

The thermodynamic potential Vt for massless quarks is given by 

Vl = -kT ^ In [1 + exp p{fXg - p)] . 

Hence 

Q kTd ( [°° ti^i \ 

= ^T7~^r^/ P dpln[l+exp/5(/i5-p) . 
Performing the last integration by parts and rearranging the terms, one finally gets 



^2 = ^ - -^{f^t 1^2 - (/i, - A)=' In [1 + exp(/3A)]}. (14) 

The quantal pressure P3 is given by 

P3 = -{dE,/dVa)Ns = Ns/SR",. (15) 

It is easy to verify that if A is set equal to zero, i.e. if no pairing is allowed to take place, N3, 
E3, P3 vanish and N2, E2, S2, P2 become equal to Ni, Ei, Si, Pi, respectively. Contribution 
of antiquarks to the above thermodynamic quantities is smaller by many orders of magnitude. 
Although we have included it in our calculations its effect is negligible. 
The droplet pressure P^ is given by 

Pd = P2 + P3 + P2 + P3 + Pgluons, (16) 

where P2 and P3 are as in Eqs. (14)-(15), P2 and P3 are the corresponding pressures due to 
antiquarks and Pgiuons is the pressure due to gluons in the droplet. 

Consider once again the spherical region of radius Rd- Initially, i.e. before the pairing 
interaction is switched on, there are A'^i quarks, Ni antiquarks and a certain number of 
gluons in it. In the final state, i.e. when the pairing interaction is on, there are N2 quarks, 
{N3/2) quark pairs, N2 antiquarks, {N3/2) antiquark pairs and the same number of gluons 
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as before. The activation energy (A) is defined as the difference between the Helmholtz free 
energies {F) in the final and initial states. Since the gluonic contribution cancels out, one 
gets 

A = FiN2) + FiN3/2) - FiN,) + ^(iVa) + F(iV3/2) - F(iVi) 

= E2- TS2 + E3-Ei + TSi + E2- TS2 + E3-E1 + TSi. (17) 

In the limit A ^ 0, this A vanishes. 

The reader may wonder why surface tension of the droplet makes no appearance in our 
formalism. In the classical theory of nucleation, there are two competing contributions to 
the energy of a droplet: droplet formation reduces bulk free energy but only at the cost of 
surface energy. The former contribution wins if the droplet radius exceeds a certain critical 
value, while the surface tension ensures that a droplet with a smaller radius shrinks and 
disappears. In the formalism presented here, energy of the droplet is calculated quantum 
mechanically. It is clear from Eq. (9) that small droplets have large quantal energy and 
tend to dissolve back into the QGP gas. This is similar to the effect of the surface tension 
in the classical theory. Moreover, the thermal wavelength (A) of the quarks is so large that 
a classification of particles as those localized on the surface and those in the interior breaks 
down. In short, we did not neglect the surface energy; it is included in the total free energy 
of the quantal droplet. 

III. RESULTS AND DISCUSSION 

The precise boundaries of the various phases in Fig. 1 are not known. Since the baryon 
number density (ub) of the normal nuclear matter is ~ 0.17 fm~^, we have performed numer- 
ical calculations for ub = 0.5 to 3 fm~^ and temperature T = to 100 MeV. The resulting 
quark chemical potential fig varies over ~ 250 — 600 MeV and the baryon chemical potential[| 



^This definition is applicable only for a gas of quarks, and not for the normal nuclear matter. 



fiB = Sf^q varies over ~ 750 — 1800 MeV. This fully covers the densities likely to be seen in 
the neutron-star cores. We have considered four different values of the superconducting gap 
parameter A, namely, 1, 10, 50 and 100 MeV. 

Since we are interested in a gas of weakly interacting quarks, antiquarks and gluons, we 
take as to be zero. Small values of ag do not affect the results in Eqs. (l)-(3) significantly. 
We first choose a set of values for ns and T, and solve Eq. (1) numerically to determine fiq. 
This fig is then used to evaluate Pgas as given by Eq. (3). Also, once fig is determined, Uq 
and rig can be evaluated separately by using 

(Above integrations cannot be done analytically, and so it is not possible to get closed 
expressions like Eq. (1), for Ug and riq.) Next, we choose a value for A, and treating the 
droplet radius Rd as a free parameter, calculate various other quantities; see Eqs. (4)- (17). 
Some of these results are shown in Figs. 3-4. 

In our formalism, as in the classical theory, the activation energy when plotted as a 
function of the radius (Fig. 3), initially increases from zero, attains a maximum, then starts 
decreasing monotonically eventually becoming negative for a sufficiently large value of the 
radius. In other words, the activation energy exhibits a barrier. The barrier originates as 
follows. Consider Eq. (17). As already stated, the antiquark terms are extremely small. 
The terms TSi and TS2 are also negligible and in fact they tend to cancel each other. That 
leaves three terms, namely Ei, E2 and E^. For small radii, E2 and the quantal energy 
term E^ dominate pushing the activation energy curve upwards. For large radii, the Ei 
term dominates pulling the curve to the negative domain. As expected, the barrier becomes 
wider if T is raised or if ub is lowered. The critical radii can be read off from Fig. 3. 
Droplets with larger radii readily grow because that tends to reduce the activation energy. 
This in turn easily brings about the transition from QGP to the superconducting phase. 

Figure 4 shows the ratio of the droplet pressure to the pressure of the surrounding gas. 



for A = 100 MeV. This ratio is significantly larger than unity for small droplets, which 
will tend to push the radii to larger values, thereby helping to bring about the QGP to 
superconducting phase transition. In the limit A ^ 0, this ratio tends to unity and is thus 
incapable of pushing the droplet radii to larger values. 

It is interesting to know the properties of the droplet when the pressure inside the droplet 
(Eq. 16) equals the pressure of the surrounding QGP gas (Eq. 3). These results are shown 
in Figs. 5-8. The condition P^ = Pgas fixes Rd- We solve the equation Pd = Pgas numerically 
to determine Rd, for various values of Hb, T and A. ns was varied from 0.5 to 3 fm~^ in 
the steps of 0.5 fm~^. Figure 5 shows the results for Rd- As T and/or ns are increased, Rd 
decreases. However, Rd is found to be rather insensitive to the variation in hb- Although we 
have varied T between and 100 MeV, temperatures relevant for neutron stars are close to 
zero. Because of the small value of Rd, the energy of even the lowest excited state inside the 
droplet turns out to be several hundreds of MeV which is large compared to the temperature, 
effectively forcing all Cooper pairs (bosons) to occupy the ground state. Note also that at low 
temperatures, the droplet radius Rd is larger than the Pippard coherence length ^ ~ l/(7rA). 

Figure 6 shows the activation energy A. First see the upper half of Fig. 6. For a fixed 
Hb, as T increases, A changes its sign from negative to positive. This is expected because, 
at higher T, droplet formation is energetically unfavourable. Secondly, as hb increases (see 
curves labelled ato f), A remains negative up to higher and higher T. This is also expected 
because at higher ub, the transition to the superconducting phase is facilitated. Now see 
the lower half of Fig. 6. It shows the A = contour, in the T — hb plane. This is also the 
phase transition curve separating the QGP phase at high T and the superconducting phase 
at low T: Recall A = Fsc — Fqgp where F are the free energies of the two phases. Both 
Fsc and Fqgp decrease as T increases and cross each other at some temperature, say T^ 
(Fig. 7). Since the system follows the path of minimum F , there is a kink in F{T). In other 
words, the slope dF/dT = —S has a discontinuity, indicating a first-order phase transition 
at Tc- Thus the A = contour serves to demarcate the region of the superconducting phase. 
Note that as np increases, the A < region in Fig. 6 survives up to higher and higher T. 



We now present an estimate of the latent heat (L) associated with the first-order phase 
transition considered in this paper. The transition temperature Tc for any given ub can be 
determined from Fig. 6. The latent heat is then estimated as L = Tc x (Si{Tc) — 5*2 (Tc)). 
For A = 100 MeV, this turns out to be ~ 150 to ~ 250 MeV and it increases with ub- In 
the limit A — i> 0, the two free-energy curves in Fig. 7 coincide, A vanishes, 5*1 = S2, there 
is no discontinuity in entropy, and the latent heat is zero. There is only one phase namely 
QGP, and there is no question of any phase transition. 

Finally, Fig. 8 shows 3-dimensional surface plots of the (normalized) activation energy 
A' defined as A' = {Fsc — Fqgp) / {Fsc + Fqgp) as a function of hb and T, for A = 1 and 
100 MeV. It is clear that in the limit A ^ 0, the activation energy surface is flat and there is 
no energy advantage for the QGP to superconducting phase transition. Indeed there would 
be only one phase, namely QGP, which would extend up to the T = axis. On the other 
hand, for A = 100 MeV, the surface is such that the system would like to "roll down" the 
slope thereby bringing about the phase transition. 

It is interesting to compare the activation energy results for the QCD phase transition, 
obtained here with those for the Bose-Einstein condensation in a gas of trapped atomic 
hydrogen, presented in |Q. In |^ two gas densities were considered: 10^^ cm~^ and 7 x 10^^ 
cm^^ both at the same temperature 50 ^K. We have repeated their calculations for T = 
45 — 55 ^K. Their activation energy exhibits a barrier, and they denote the maximum value 
of the activation energy by Aq. Consider the ratio A^/kT at the peak position and the 
corresponding number A'"o of particles in the droplet. In these two quantities are highly 
sensitive to the changes in the gas density and temperature.0 For instance, for the above 
tiny change in the gas density, the ratio Ao/kT changes from 4.8 to 104.1, and Nq changes 
from 18 to 3862. Similar sensitivity was seen when T was changed by a few fiK. In contrast. 



^There is a misprint in j^: in Fig. 1, the y-axis label should read A{10 ^^) a.u. For the inset, it 
is ^(10"'^) a.u. 
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the results presented here are relatively stable with respect to changes in ns and T. 

In conclusion, we have studied hitherto unexplored first-order phase transition from QGP 
to a colour-superconducting phase, and have presented a simple model of the transition 
with very few assumptions. It provides numerical support to the prevalent ideas of phase 
structure of QCD in the low-temperature, high-baryon-number-density regime. We have 
studied properties of the droplets of the superconducting phase, formed at various points in 
the T — hb plane, for several values of the gap parameter. We have been able to estimate 
the extent of the superconducting phase in this plane and the latent heat of the transition. 
At low temperatures, droplet radii are found to be larger than the Pippard coherence length 
which measures the spatial extent of the qq Cooper pairs. 

RSB acknowledges the hospitality of the McMaster University where this work was ini- 
tiated. We thank Mark Alford and Mustansir Barma for their comments on this work. This 
research was partially supported by a grant from NSERC (Canada). 
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Appendix 

If f{p) is tlie distribution function of a gas in equilibrium, the number density n is given 
by ^ = / f {p)dP'p / {^-kTiY . For a nonrelativistic classical gas at temperature T, one has 



A^-^] 



f = exp 

where f3 = 1/kT and n is the chemical potential. Hence 

n = exp{f3fi){27rh^ /mkT)~^/^ = exp(/3/i)A-^ 

where A is the thermal wave length. 

For an extreme relativistic gas of massless classical particles, 

/ = exp [-/3{pc - /i)] , 

and hence 

n = exp(/3/i)(c/i7r2/7A;T)-3. 

We define the thermal wave length in this case as A = chn^'^ /kT, so that the fugacity 
exp(/9/i) once again equals n\^. 
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FIGURE CAPTIONS 

FIG. 1. Schematic phase diagrams of QCD: T is the temperature and fis is the baryon 
chemical potential. E is the critical point at which the line of first-order phase transitions 
ends. 2SC and CFL (Colour-Flavour-Locked) are 2-fiavour-like and 3-flavour-like colour 
superconducting phases, respectively. The dashed line denotes the critical temperature 
at which quark-quark pairing vanishes. The thick arrow indicates the first-order phase 
transition considered in this paper. For the normal nuclear matter at T = 0, fiB should be 
defined as mjv + h'^kp/2mN where kp ~ 1.4 fm"-*^ |^ and mj^ is the nucleon mass, which 
gives fiB — 980 MeV. (Figure adapted from [l|].) 

FIG. 2. Feynman diagrams contributing to the equation of state of quark-gluon plasma 
in order a^- Curly lines: gluons, dashed lines: ghost particles, solid lines: quarks. 

FIG. 3. Activation energy A vs droplet radius Rd for various gap values A. Baryon 
number density ub (fm~^) and temperature T (MeV) are: (a) 0.5 and 5 (b) 0.5 and 100 (c) 
3 and 5 (d) 3 and 100, respectively. 

FIG. 4. Ratio of the pressure of the droplet to the pressure of the surrounding QGP gas 
vs droplet radius Rd, at gap A = 100 MeV. Curve labels are as in Fig. 3. 

FIG. 5. Droplet radius Rd vs temperature T for various gap values A. Results are 
insensitive to ub, so not all curves are labelled. 

FIG. 6. Upper panel: Activation energy A vs temperature T. Baryon number densities 
Ub (fm~^) are (a) 0.5 (b) 1 (c) 1.5 (d) 2 (e) 2.5 (f) 3. Lower panel: The A = contour or 
the phase transition line in the T — ub plane. SC: superconducting phase. 

FIG. 7. Free energies Fqgp and Fsc in the QGP and superconducting phases, respec- 
tively, vs temperature T. 

FIG. 8. Three-dimensional surface plots of the activation energy A' for A = 1 MeV 
(upper part) and A = 100 MeV (lower part), as a function of ub (fm^'^) and T (MeV). 



14 



T 



150 
MeV 



mg 


» mu,d ^ 




E 


QGP 




Hadronic 







Matter 


Y''^2SC 


^CFL 




Hadronic 
Matter 



CFL 



Nuclear 
Matter 



1 GeV 



^^B 



FIG. 1. 




+ 




+ tnnre~s~s~s-i + 




FIG. 2. 



■200 - 



200 - 



100 - 



> 

0) 



A=10MeV- 




A=100MeV- 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 



FIG. 3. 






Q_ 




0.5 



1 1.5 



2 2.5 



8 
6 



T3 4 



2 h 




"I 1 1 1 1 1 1 1 r 



A=10MeV 



J I I I I I I I L 



A = 50 MeV 




■a 
CC 



2 - 



1 - 







"1 1 1 1 1 1 1 1 r 



rig = 0.5 fm 



A= 100 MeV 




rig = 3 fm 



J I I I I I I I L 



10 20 30 40 50 60 70 80 90 100 

T (MeV) 



FIG. 5. 



300 



200 - 



> 

CD 




■200 - 



40 50 60 70 80 90 100 

T (MeV) 



> 

0) 



90 


1 1 


1 1 1 
QGP ^^^^-^--"'^ 


80 
70 


A>0 ^ 


^""sc 


60 


^^ 


- 


50 


/ A<0 


- 


40 

on 


1 1 


A= 100 MeV - 

1 1 1 



0.5 1 1.5 



2.5 



FIG. 6. 



> 





L_ 
(D 

C 

LU 

CD 

CD 

v_ 



10^ - 



10^ - 



A=100MeV 
rig = 0.5 fm' 




10' 







20 



40 60 

T (MeV) 

FIG. 7. 



80 



100 



0, 
0, 
0, 

■0, 
■0, 
■0, 
■0, 



0, 
0, 
0, 

■0, 
■0, 
■0, 
■0, 



3 
2 
1 

1 
2 
3 
4 



3 
2 
1 

1 
2 
3 
4 




